Current noise of a resonant tunnel junction coupled to a nanomechanical oscillator 
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We present a theoretical study of current noise of a resonant tunnel junction coupled to a nanome- 
chanical oscillator within the non-equilibrium Green's function technique. An arbitrary voltage is 
applied to the tunnel junction and electrons in the leads are considered to be at zero temperature. 
The properties of the phonon distribution of the nanomechanical oscillator strongly coupled to the 
electrons on the dot are investigated using a non-perturbative approach. An analytical calculations 
and numerical results for the current-voltage, shot noise and the corresponding Fano factor as a 
function of applied bias show significant features of the nanomechanical oscillator coupling dynam- 
ics. This will provide useful insight for the design of experiments aimed at studying the quantum 
behavior of an oscillator. 

PACS numbers: 73.23.Hk,85.85.+j 



I. INTRODUCTION 

In recent years, there has been a great interest in the 
measurement of current noise near the quantum limit of 
displacement sensitivity using single electron transistors 
(SET) and nanomechanical oscillators « . Current noise 
is the non-equilibrium fluctuation which is caused by the 
discreteness of charged carriers. When the size of an 
electromechanical system reaches the nanometer scale, 
the current noise problem becomes a very interesting as- 
pect of the NEMS£"i2i based devices. In NEMS devices 
it is a fundamental physical signature, which can provide 
us additional information regarding quantum transport, 
in addition to differential conductance and current volt- 
age characteristics. The electrons traveling through the 
device become correlated in the same channel and the 
same probe, as well as in different channels and different 
probes. For a mesoscopic system, the electrons are cor- 
related due to coherent transport and they are governed 
by the Fermi distribution and Coulomb blockade. Noise 
is caused by the randomness of electron scattering and it 
has been extensively studied in different types of meso- 
scopic structures. For uncorrelated electrons travelling 
through a macroscopic conductor, shot noise 5(0) at zero 
frequency is given by the Poisson value S(0) = 2el, where 
I is the net average current flowing through the system. 
However, for a mesoscopic system, electrical current cor- 
relations are dominated by the Coulomb blockade and a 
Fermi distribution can lead to the deviation of shot noise 
from a Poissonian (F=l) form. An important parameter 
for describing this type of shot noise is the Fano factor 
F= S(0)/2eI, using which the shot noise can be classified 
as sub- Poissonian (F < 1) or super-Poissonian (F > 1). 
Similarly, phonon assisted features have been observed in 
molecular systems^— with strong interactions, which is 
beyond the scope of perturbation theory. 

In general, there are two different theoretical formu- 
lations that can be used to study the quantum trans- 
port in nanoscopic systems under voltage bias. Firstly, a 
generalized quantum master equation approach^— and 
secondly, the non-equilibrium Green's function formula- 



tionSIr— . The former leads to a simple rate equation, 
where the coupling between the dot and the leads is con- 
sidered as a weak perturbation and the electron-phonon 
interaction is also considered very weak. In the lat- 
ter case one can consider strong leads to system and 
electron-phonon coupling. The non-equilibrium Green's 
function technique is able to deal with a very broad va- 
riety of physical situations related to quantum transport 
at molecular— ~— levels. It can deal with strong non- 
equilibrium situations and very small to very large ap- 
plied bias. In most of the theoretical work on NEMS 
devices, the mechanical degree of freedom has been de- 
scribed classically/semiclassicall y 15 ' 16 using the master 
equation approach or quantum mechanically^£ ~ 20 i 34 i 35 
using the perturbation approach for the electron-phonon 
interactions. Furthermore, Shot noise for electron trans- 
port in molecular devices was investigated within a scat- 
tering theory approac h 29 i 30 . Phonon effects in the noise 
spectrum were studied recently 3 -! - — , in connection with 
NEMS and both single electronic levels in quantum dots 
and single molecules. In Armour's work 15 , the me- 
chanical part was also treated classically, including the 
damped oscillator and assuming a weak electron tun- 
neling process. This approach is based on a perturba- 
tion, weak coupling and large applied bias approxima- 
tions, whereas the Keldysh non-equilibrium Green's func- 
tion formulation can treat the system-leads and electron- 
phonon coupling with strong interactions for both small 
and large applied bias voltage. Moreover, these the- 
ories fail to explain the low bias regime. The trans- 
port properties have been described and discussed semi- 
classically/classically but need a complete quantum me- 
chanical description. A theory beyond these cases is re- 
quired to further refine experiments to investigate the 
quantum transport properties of NEMS devices. 

In the present work, we employ the non-equilibrium 
Green's function method to discuss the current-voltage 
and shot noise properties of a NEMS device. This is a 
fully quantum mechanical formulation whose basic ap- 
proximations are very transparent, as the technique has 
already been used to study transport properties in a wide 



range of systems. The main differences between exist- 
ing work and ours are: in most of the existing literature 
a very large chemical potential difference is considered 
while we are able to include a range from very small to 
very large. In our calculation the inclusion of the oscilla- 
tor is not perturbative as the STS experiments^^— are 
beyond the range of perturbation theory. Hence, an ap- 
proach is required beyond the quantum master equation 
or linear response. Hence, our work provides an exact an- 
alytical solution to the current-voltage, shot noise, cou- 
pling of leads with the system, very small to very large 
chemical potential difference and includes both the right 
and left Fermi level response regimes. For simplicity, we 
used the wide-band approximation, where the coupling 
between the leads and the dot is taken to be indepen- 
dent of energy. This provides a way to perform transient 
transport calculations from first principles while retain- 
ing the essential physics of the electronic structure of the 
dot and the leads. Another advantage of this method 
is that it treats the infinitely extended reservoirs in an 
exact way in the present system, which may give a bet- 
ter understanding of the essential features of NEMS in a 
more appropriate quantum mechanical picture. 
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where eo is the single energy level of electrons on the dot 
with Cq,Co the corresponding creation and annihilation 
operators. The parameter r\ = \lj\[2 physically repre- 
sents an effective electric field in the capacitor formed 
by the oscillator and the dot electrons , which we shall 
call coupling strength between the oscillator and the elec- 
trons on the dot given as A^ = eEl, where e is the charge 
of electron, E is the strength of the electric field and 

I = \ — is the zero point amplitude of the oscillator 

with mass fi. Here we assume that the energy of the res- 
onant level depends linearly on the oscillator coordinate. 
The frequency of the nanomechanical oscillator is u> and 
b\ b are the raising and lowering operator of the phonons 
given as 



II. MODEL CALCULATIONS 

We consider a single quantum dot connected to 
two identical metallic leads via tunnel junctions (SET 
device) 15 ' 19 i 35 . A single nanomechanical oscillator is cou- 
pled to the electrons on the dot and the applied gate 
voltage is used to tune the single level of the dot. This 
model represents a close analog of an electromechanical 
displacement detector and contains the essential features 
of the recent experiments performed with a resonatori 2 ^, 
or a vibrating beam of crystal or cantilever coupled to 
the SET device 3-6 . In most of these experiments, the 
motion of a resonator (nanomechanical oscillator) may be 
detected by capacitively coupled electrodes placed on the 
resonator and biasing the electrode at a constant voltage. 
The dot has a small capacitance due to its small diam- 
eter (nm) and thus has charging energy. This charging 
energy exceeds the thermal energy in these experiments. 
For this reason we consider that only one excess electron 
may occupy the device. To ensure this condition we are 
working at zero temperature. 

In the present simple model, both the electronic and 
the mechanical degrees of freedom have been treated 
quantum mechanically, where the nanomechanical oscil- 
lator is represented in terms of creation and annihila- 
tion operators. We neglect the spin degree of freedom 
and electron-electron interaction effects and consider the 
simplest possible model system. We also neglect the ef- 
fects of finite electron temperature of the lead reservoirs 
and damping of the oscillator. Our model consists of the 
individual entities such as the single quantum dot with 
perpendicular oscillator on it and the left and right leads 
in their ground states at zero temperature. The Hamil- 
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The remaining elements of the Hamiltonian are 
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where N is the total number of states in the lead, V a 
is the hopping between the dot and the leads a — L,R, 
j represents the channels in one of the leads. For the 
second lead the Hamiltonian can be written in the same 
way. 

The total Hamiltonian of the system is thus H = Ho + 
AH a . We write the eigenvalues and the eigenfunctions 
of -ffdot-ph as 



= e + M™ + 1) -A 
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A n exp[- l ^-]H n (lK), (9) 
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for the occupied, Xq ^ and unoccupied, xq — 0, dot 
respectively, where A n = l/\j\/n2™nH, A = A 2 //zw 2 . 
xq = X/fiu> 2 is the shift of the oscillator due to the cou- 
pling to the electrons on the dot and H n (lK) are the 



usual Hermite polynomials. Here we have used the fact 
that the harmonic oscillator eigenfunctions have the same 
form in both real and Fourier space. 

In order to transform between the representations for 
the occupied and unoccupied dot we require the matrix 
with elements 



The particle current I a into the interacting region from 
the lead is related to the expectation value of the current 
operator I a = Tr(pJ a )^r—, where 

Ja = -r y^^{Vo, a clc a j - c' a;j coVo, a } , (14) 
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which may be simplified as 
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(12) 
^f- and i™ "(x) are the associ- 
ated Laguerre polynomials. Note that the integrand is 
symmetric in to and n but the integral is only valid for 
n < to. Clearly the result for n > to is obtained by 
exchanging to and n in equation (|1 2|) to obtain 
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The position of the resonant level with respect to the 
chemical potential in the leads is thus affected by the dis- 
placement of the nanomechanical oscillator, which in turn 
affects the transport properties of the junction through 
the device. 



p = —iG < and p is the density matrix written in terms 
of the lesser Green's function. Eq. (TH| can be written in 
terms of lesser Green's function as 



la® = j:{G^ a (t,t)V afi (t) - V * a (t)G< (t,t)}, (15) 
where we have the following relations 
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and 



G< (M) - jdt'{g r a Jt,t')V afi (t')G< (t',t) 

+9la(t,t')V a , (t')Gl (t',t)}, (17) 

where ga,S v-^') refers to the unperturbed states of 
the leads. 
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Using the fact that 

£^ (<) (i',*) = KL Wi a) ' (<) (^)Vko(i), 
we can simplify the above equation as 
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( 1Q ) 
In the stationary regime all of the Green's functions depend on time order (t — £'), yielding the Fourier/Laplace 

transform of equation (|19j) as 
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In order to calculate the analytical results and to dis- cuss the numerical quantum dynamics of the nanome- 



chanical system, our focus is firstly to derive an analyti- 
cal relation for the retarded self-energy The self-energy 
represents the contribution to the dot energy, due to in- 
teractions between the dot and the leads. In obtaining 
these results we use the wide-band approximation where 
the retarded self-energy of the dot due to each lead is 
considered to be energy independent and is given by 



K .,n jE) = AH* a g r a>a (E - (no + \)hw) AH a , (21) 
where off- diagonal element of matrix, E r / (E). are 
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zero and g r a a (E — (no + §)&*;) is the uncoupled Green's 
function in the leads as 
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Nn a de a , j stands for every channel in 



each lead and n a is the density of states in lead a. With 
the help of equation (1221 . the retarded self-energy may 
be written as 
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which is now independent of E, T a — 2tt \Vo, a I n a , a rep- 
resenting the L or R leads and the retarded self energy is 
now independent of the oscillator's index (no, no). Hence, 
it can be written as K n J E ) = (K .n ,a(E))* = 

2 ' 

We solve Dyson's equation using Yidot-iead as a per- 
turbation. For the more general systems we aim to treat 
in the future, this is a reasonable small parameter. In 
the present case, however, we can find an exact solution. 
The retarded and advanced Green's functions on the dot, 
with the phonon states in the representation of the un- 
occupied dot, may be written as 



G r Jk ) (E) = y / ^n, m g^ a) (E)^ not 



where gm (E) is the retarded (advanced) Green's func- 
tion on the occupied dot, 

g r J a) (E) =[E-e - (m + i)fiw + A ± irf ' , (25) 

assuming that Tl = Tr = T. 

The lesser Green's function in the presence of the 
nanomechanical oscillator including the dot and the leads 

is written as 



G^AE) = £ Gl no (E)X< on ,(E)G^ n ,(E), (26) 

n ,n' 

with S < , (E) being the the lesser self energy which is 
given as 
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where the off-diagonal element of matrix, S 
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arc 1 



zero and the diagonal (n = no) element of the lesser 
self-energy may be written as 

= ir Q / de a f a (e a ) B no 5(E -e a - (n + \)fw) 
= iT a f a (e)B no> (28) 

where e = E — no + (h)fvuj, f a (e a ) and f a {s) are the 
Fermi distribution functions of the left (a — L) and right 
(a = R) leads, which have different chemical potentials 
under a voltage bias and B no is the Boltzmann factor for 
the oscillator state. Index no determines the statistical 
occupation probability of the phonon state |no) at finite 
temperature and therefore the accessibility of particular 
conduction channels is determined by a weight factor of 
the Boltzmann distribution function. 

Following equation (I2U1) . the formula for the current 
through each of the leads is written in terms of oscillator 
indices as 
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With the help of equation (|29[) , the net current through 
the dot and the leads with the oscillator on the dot is 
written as 
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The resulting expression for the net current is 
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(31) 



which is derived from equation (|30ft using the same damp- 
ing factor for each lead (Yl =Tr = T). 

For the present case of zero temperature the lesser self- 
energy may be recast in terms of the Heaviside step func- 
tion 9(x) as 



n ,n a ,a 



\T a 9 (e Fa - e) 5. 
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where e Fa is the Fermi energy on lead a and the Kro- 
necker delta, 5 n0r o, signifies that the nanomechanical os- 
cillator is initially in its ground state, n® = 0. 

Inserting this equation in the above equation, the re- 
sult is 



After performing the integral in the above expression, the 
final result is written as 
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where 
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T (E) = -(G r 0fi (E) - G a 0t0 (E)) = T[G r 0t0 (E)TGl (E)] 

Using the expression for the retarded and advanced 
Green's function, the expression for the net current be- 
comes 
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The model represent the interplay between two physical 
time scales of the system, the oscillator frequency and the 
tunneling rate. This model also shows the very interest- 
ing interplay between two physical length scales of the 
system, zero point amplitude and zero point displace- 
ment, which actually affected by the weak and strong 
coupling dynamics. 



The zero frequency shot-noise has been derived and 

applied successfully in many examples of the transport 

(34) dynamics of nanoscopic system s 17 ' 24 i 26 . This is given as 
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Using the values of transmission coefficients at zero temperature, the above expression can be simplified as 
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After integrating the above expression, we arrive at the final result 

where S n (Q) and 5 n>m (0) are defined as 
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The corresponding Fano factor can be calculated from the 
zero frequency noise and the net current flowing through 
the system, which is defined as 



III. DISCUSSION OF RESULTS 

The current-voltage (I-V) characteristics, shot noise 
spectrum and the corresponding Fano factor of a reso- 
nant tunnel junction coupled to a nanomechanical oscilla- 
tor are shown graphically for different values of coupling 
strength, using the same parameters as: the single energy 
level of the dot eo = 0.5, the characteristic frequency of 



tan 



e F fl — eo — mhw + A 
f 



r 



the oscillator hw = 0.1, the damping factor T = Q.lhw 
and the chemical potentials < e F ^ < 1 and e F # = 0. 
These are chosen to illustrate the physics of such systems 
rather than to represent a specific implementation. The 
oscillator induced resonance effects are clearly visible in 
the numerical results. It must be noted that we have ob- 
tained these results in the regime of strong and zero or 
weak coupling of the oscillator with the electrons on the 
dot. The coupling between the leads and the dot is con- 
sidered to be symmetric and we assume that the electrons 
in the leads are at zero temperature and have constant 
density of states. With increasing coupling strength, the 
number of additional steps also increases while for zero or 
weak coupling we find only the basic resonance step. This 
confirms the effect of the coupling between the electrons 
on the dot and the single oscillator mode where higher 
energy electrons are able to drop to the dot energy by 



creation of phonons. 

Closer analytical examination of the expression for the 
current and shot noise (35 & 38) shows that the main res- 
onance steps occur when the applied voltage, 6fl is equal 
to the energy eigenvalues of the coupled dot electron and 
the nanomechanical oscillator. The main step (n = 0) 
is given by the Lorentzian form with its center at the 
eph — £o — A, known as a Breit-Wigner resonance. The 
additional steps due to the emission of phonons can be 
seen on the positive energy side with tpL — e — A + nHuj 
where w is the characteristic frequency of the nanome- 
chanical oscillator. We have shown the I-V charactcris- 
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FIG. 1: Current-voltage characteristics (dimensionless) as a 
function of applied voltage €fl (in arbitrary units), coupling 
strength r\ — 0.1fia;(dotted line) and lftio (solid line). Gate 
voltage en = 0.5, oscillator energy fhu — 0.1 and self-energy 

r = o.ifiw. 

tics of the NEMS device against applied bias for differ- 
ent values of the coupling strength in Fig. [1] The main 
resonance step is the elastic or zero phonon transition. 
The amplitude of the additional steps is much smaller 
than the basic resonance step. The electrons that tunnel 
onto the dot can only excite the oscillator mode as at 
zero temperature there are no phonons available to be 
absorbed. Moreover, we have seen that with increasing 
coupling strength, the number and intensity of the addi- 
tional steps increases but their intensity always remains 
much smaller than the main step. The steps in the cur- 
rent characteristics vanish if the upper electrochemical 
potential is smaller than the dot energy plus the oscilla- 
tor energy. 

Next, we have shown the shot noise as a function of 
applied bias in Fig. [5J which exhibits the single step in 
the presence of weak or zero electron-oscillator coupling 
while the number of additional steps increases with in- 
creasing coupling strength. Obviously, the shape of the 
shot noise curve is similar to that of the net current, 
as shown in Fig. [TJ The only difference is associated 
with their behavior above the resonance point, where the 
noise power for the strong coupling case can exceed the 
shot noise for the zero phonon case. We also show the 
differential shot noise against applied bias for different 
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FIG. 2: Shot noise (dimensionless) as a function of applied 
voltage 6fl (in arbitrary units), with gate voltage en=0.5, os- 
cillator energy fhu — 0.1, self-energy T = O.lfej, and coupling 
strength 77 = O.lfej (dotted line) and lftui (solid line). 
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FIG. 3: Differential shot noise (dimensionless) as a function 
of applied voltage €fl (in arbitrary units) with gate voltage 
en = 0.5, oscillator energy fku — 0.2, self-energy T — 0.2fej 
and coupling strength r\ = O.lfiw (dotted line) and r\ = 0.5fej 
(solid line). 



values of the coupling strengthen Fig. [3]). In the ab- 
sence of phonons, when the transport is coherent, the 
shot noise spectrum exhibits two peaks separated by an 
antiresonance and located symmetrically around the po- 
sition where the current step associated with the single 
level of the dot is located. The origin of such an antires- 
onance in the noise spectrum is associated with the fact 
that no noise is generated when the transmission via the 
dot state is perfect, T = 1 or zero. Due to the presence of 
the nanomechanical oscillator, the main resonance peaks 
are shifted by A. Meanwhile, the sharp peaks are now 
accompanied by a set of additional peaks. We note that 
the separation between the differential shot noise peaks is 
set by the frequency of the oscillator. This phenomenon 
can be explained on the basis of phonon emission during 
electron-phonon coupling, indicating that a new channel 
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has opened and contributes to the resonant tunneling 
process. 



1.0 
0.9 




i V 5 


: 


0.8 


: 




■ 


0.7 


]/" : 


0.6 


: 




: : 


<i 


: 


: 




; .-" 


0.4 





■* 



0.0 



0.2 0.4 0.6 

Applied Voltage (Arb. Units) 



0.8 



1.0 



FIG. 4: Fano factor as a function of applied voltage efl (in 
arbitrary units) with gate voltage eo=0.5, oscillator energy 
Ikd — 0.1, self-energy T = O.lfej and coupling strength 77 = 
0.1/zw (dotted line) and r\ = lhcu (solid line). 

Finally, we have shown the Fano factor against ap- 
plied bias for different values of the coupling strength. 
Information about the statistical properties of the elec- 
trons is included in the Fano factor, which is plotted in 
Fig. 21 Since in our model all the interactions between 
the current carriers are neglected, such electron corre- 
lations are associated only with the Coulomb blockade. 
This principle is related to the fact that one electron 
feels the presence of the others, since it cannot occupy 



the state on the dot already occupied by the electron. 
The crossover in the shot noise power from Poissonian 
(F =1) to sub-Poissonian (F<1) is always observed after 
the first step in the current voltage dependence. This 
implies that electrons tunnel in a correlated way in the 
NEMS device. The most important result is the signifi- 
cant enhancement of the Fano factor due to the phonon 
effects, observed for cfl > £o — A + nhu>, where the 
multi-channel process reduces electron correlations com- 
pared with the single channel case. Moreover, the shift 
A can also be easily recognized in Fig. |4] 



IV. SUMMARY 

In this work, we analyzed the current noise character- 
istics of a resonant tunnel junction coupled to a nanome- 
chanical oscillator by using the non-equilibrium Green's 
function approach without treating the electron-phonon 
coupling as a perturbation. We have derived an analyt- 
ical expressions for the net current flowing through the 
system and for the shot noise. This enables us to see the 
effects of the coupling of the electrons to the oscillator 
on the dot and the tunneling rate of electrons. We show 
the numerical results for very weak and strong coupling 
strength. We have found additional steps or peaks due to 
coupling of single phonon mode which are absent for very 
weak or no coupling strength. We also discuss the corre- 
sponding Fano factor as a function of applied bias which 
shows thermal or poissonian behavior to non-thermal or 
sub-poissonian behavior. 
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